/* Benegal and Scruggs misinformation correction study on renewables
Replication file */

clear
** set wd
use RRdata
set more off


// FIGURE 2

eststo clear
eststo rel1: reg reliablescaled b2.treatment, vce(hc2)
eststo rel3: reg reliablescaled b2.treatment if Republican==1, vce(hc2)
eststo rel4: reg reliablescaled b2.treatment if Democrat==1, vce(hc2)

set scheme plottig
coefplot (rel1, label (All respondents) msymbol(X)) (rel3, label (Republicans) msymbol(S) mcolor(red) ciopt(lc(red))) (rel4, label (Democrats) msymbol(T) ciopt(lc(blue)) ///
mcolor(blue)), drop(_cons) xline(0) xtitle("Difference vs. uncorrected misinfo treatment.")


# delimit ;
esttab using reliability.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) r2(%9.3f) 
	title("Reliability of solar and wind energy" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("OLS with HC2 robust standard errors")
	compress replace;
 #delimit cr
 
set scheme plottig
coefplot (rel1, label (All respondents) msymbol(X)) (rel3, label (Republicans) msymbol(S)) (rel4, label (Democrats) msymbol(T)), ///
drop(_cons) xline(0) xtitle("Difference vs. uncorrected misinfo treatment.")

// FIGURE 3

// how much from:: solar (3), wind (4)
eststo clear
eststo solar1: reg solarscaled b2.treatment, vce(hc2)
eststo wind1: reg windscaled b2.treatment, vce(hc2)

eststo solar2: reg solarscaled b2.treatment if Republican==1, vce(hc2)
eststo wind2: reg windscaled b2.treatment if Republican==1, vce(hc2)

eststo solar3: reg solarscaled b2.treatment if Democrat==1, vce(hc2)
eststo wind3: reg windscaled b2.treatment if Democrat==1, vce(hc2)

coefplot (solar1, label (All respondents) msymbol(X)) (solar2, label (Republicans) msymbol(S) mcolor(red) ciopt(lc(red))) (solar3, label (Democrats) ///
mcolor(blue) ciopt(lc(blue)) msymbol(T)), drop(_cons) xline(0) title("Solar") xtitle("Difference vs. uncorrected misinfo") legend(rows(1))
graph save solar, replace

coefplot (wind1, label (All respondents) msymbol(X)) (wind2, label (Republicans) msymbol(S) mcolor(red) ciopt(lc(red))) (wind3, label (Democrats) ///
mcolor(blue) ciopt(lc(blue)) msymbol(T)), drop(_cons) xline(0) title("Wind") xtitle("Difference vs. uncorrected misinfo")
graph save wind, replace

grc1leg solar.gph wind.gph, col(2) ycommon xcommon graphregion(margin(none)) legendfrom(solar.gph) ring(1)

eststo clear
eststo solar1: reg solarscaled b2.treatment, vce(hc2)
eststo solar2: reg solarscaled b2.treatment if Republican==1, vce(hc2)
eststo solar3: reg solarscaled b2.treatment if Democrat==1, vce(hc2)

# delimit ;
esttab using solarols.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) r2(%9.3f) 
	title("Preferred usage" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("OLS with HC2 robust standard errors")
	compress replace;
 #delimit cr
 
eststo clear
eststo solar1: oprobit solarscaled b2.treatment
eststo solar2: oprobit solarscaled b2.treatment if Republican==1
eststo solar3: oprobit solarscaled b2.treatment if Democrat==1

# delimit ;
esttab using solarlogit.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) pr2(%9.3f) 
	title("Preferred usage" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("Ordinal logit, odds ratios")
	compress replace;
 #delimit cr
 

eststo clear
eststo wind1: reg windscaled b2.treatment, vce(hc2)
eststo wind2: reg windscaled b2.treatment if Republican==1, vce(hc2)
eststo wind3: reg windscaled b2.treatment if Democrat==1, vce(hc2)

# delimit ;
esttab using windols.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) r2(%9.3f) 
	title("Preferred usage" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("OLS with HC2 robust standard errors")
	compress replace;
 #delimit cr
 
eststo clear
eststo wind1: oprobit windscaled b2.treatment
eststo wind2: oprobit windscaled b2.treatment if Republican==1
eststo wind3: oprobit windscaled b2.treatment if Democrat==1

# delimit ;
esttab using windlogit.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) pr2(%9.3f) 
	title("Preferred usage" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("Ordinal logit, odds ratios")
	compress replace;
 #delimit cr
 

eststo clear
eststo taxes1: reg supporttaxes b2.treatment, vce(hc2)
eststo mandate1: reg REminimum b2.treatment, vce(hc2)

eststo taxes2: reg supporttaxes b2.treatment if Republican==1, vce(hc2)
eststo mandate2: reg REminimum b2.treatment if Republican==1, vce(hc2)

eststo taxes3: reg supporttaxes b2.treatment if Democrat==1, vce(hc2)
eststo mandate3: reg REminimum b2.treatment if Democrat==1, vce(hc2)


set scheme plottig

coefplot (taxes1, label (All respondents) msymbol(X)) (taxes2, label (Republicans) msymbol(S) mcolor(red) ciopt(lc(red))) (taxes3, label (Democrats) ///
msymbol(T) mcolor(blue) ciopt(lc(blue))), drop(_cons) xline(0) title("Increase taxes on fossil fuel usage") xtitle("Difference vs. uncorrected misinfo") legend(row(1))
graph save taxes, replace

coefplot (mandate1, label (All respondents) msymbol(X)) (mandate2, label (Republicans) msymbol(S) mcolor(red) ciopt(lc(red))) (mandate3, label (Democrats) ///
msymbol(T) mcolor(blue) ciopt(lc(blue))), drop(_cons) xline(0) title("Minimum renewable energy mandate") xtitle("Difference vs. uncorrected misinfo")
graph save mandate, replace

grc1leg taxes.gph mandate.gph, col(2) ycommon xcommon graphregion(margin(none)) legendfrom(taxes.gph) ring(1)



eststo clear
eststo taxes1: reg supporttaxes b2.treatment, vce(hc2)
eststo taxes2: reg supporttaxes b2.treatment if Republican==1, vce(hc2)
eststo taxes3: reg supporttaxes b2.treatment if Democrat==1, vce(hc2)

# delimit ;
esttab using taxols.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) r2(%9.3f) 
	title("taxes" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("OLS")
	compress replace;
 #delimit cr
 
eststo clear
eststo taxes1: probit supporttaxes b2.treatment
eststo taxes2: probit supporttaxes b2.treatment if Republican==1
eststo taxes3: probit supporttaxes b2.treatment if Democrat==1

# delimit ;
esttab using taxlogit.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) pr2(%9.3f) 
	title("taxes" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("logit")
	compress replace;
 #delimit cr
 
**

eststo clear
eststo mandate1: reg REminimum b2.treatment, vce(hc2)
eststo mandate2: reg REminimum b2.treatment if Republican==1, vce(hc2)
eststo mandate3: reg REminimum b2.treatment if Democrat==1, vce(hc2)

# delimit ;
esttab using mandateols.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) r2(%9.3f) 
	title("mandate" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("OLS")
	compress replace;
 #delimit cr
 
eststo clear
eststo mandate1: probit REminimum b2.treatment
eststo mandate2: probit REminimum b2.treatment if Republican==1
eststo mandate3: probit REminimum b2.treatment if Democrat==1


# delimit ;
esttab using RElogit.rtf, nogap se b starlevels(* .05 ** 0.01 *** .005) pr2(%9.3f) 
	title("re" )
	   label nonumbers mtitles("All respondents" "Republicans" "Democrats")
	addnote("logit")
	compress replace;
 #delimit cr
 